library(daewr)
library(mixexp)
dat <- read.csv("./data1sample.csv",header=T)
orgColnames <- colnames(dat)
colnames(dat) <- c("x1","x2","x3","x4","x5","x6","y")
head(dat)

## Quadratic Model
sqm <- lm(y~(x1+x2+x3+x4+x5+x6)^2,data=dat)
summary(sqm)
plot(sqm)
